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Abstract: This work investigates the possibihty of a long-hved stop squark in 
supersymmetric models with the neutrahno as the hghtest supersymmetric particle 
(LSP). We study the implications of meta-stable stops on the sparticle mass spectra 
and the dark matter density. We find that in order to obtain a sufficiently long stop 
lifetime so as to be observable as a stable i?-hadron at an LHC experiment, we need 
to fine tune the mass degeneracy between the stop and the LSP considerably. This 
increases the stop-neutralino co-anihilation cross section, leaving the neutralino relic 
density lower than what is expected from the WMAP results for stop masses '^1.5 
TeV/c^. However, if such scenarios are realised in nature we demonstrate that the 
long-lived stops will be produced at the LHC and that stop-based i?-hadrons with 
masses up to 1 Te V / c^ can be detected after one year of running at design luminosity. 
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1. Introduction 

The well known problems of the lack of naturalness of the Standard Model and its 
inability to explain dark matter have led to an expectation of hitherto unobserved 
particles and physics processes being manifest at TeV energy scales[l]. A number of 
candidate phenomenological extensions of the Standard Model have been proposed 
although, as yet, there is no experimental evidence to either falsify or confirm the 
correctness of any of these models. It is therefore prudent to consider a range of 
observables which could be used to discover and quantify any exotic physics signals 
seen at the TeV scale. Such observables include anomalous rates of jet and missing 
transverse energy production, and the observation of exotic stable^ massive interact- 
ing particles (SMPs), which is the subject of this paper. The detection of SMPs is 
experimentally challenging owing to their slow speeds and interactions in matter. In 
this work, it is shown how a stable stop can arise in supersymmetric scenarios [2] in 
which the neutralino is the lightest supersymmetric particle (LSP). The consistency 
of such scenarios with cosmological constraints is also studied. Stable stops would 

^In this paper, the terms stable and long-hved particle are taken to imply that the particle will 
not decay during its traversal of a typical detector at a collider unless stated otherwise. 
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manifest themselves in an experiment as bound hadronic states (termed i?-hadrons). 
Using the Pythia[3] generator, the signatures of MSSM-8 models which give rise to 
i?-hadrons at the LHC are investigated. 

Long-lived stops are predicted in many theoretical models[4, 5]. In the theory 
of electroweak baryogenesis [6, 7, 8] a long lived stop is predicted in the so-called 
coannhilation region[9, 10, 11]. The stop lifetime in these models does not typically 
exceed 10~^^ s, leaving no possibility of them traversing any substantial part of a 
detector at a high-energy collider experiment. Supersymmetry (SUSY) models with 
a gravitino LSP[12, 13, 14] can also give rise to a long-lived stop if it is the next-to- 
lightest supersymmetric particle (NLSP). Due to the smallness of the gravitational 
coupling these stops can have a macroscopic lifetime [15]. Lifetimes with similar 
values can be obtained from models where the right-handed sneutrino is the LSP[16], 
since the coupling is proportional to the neutralino Yukawa coupling. A stop NLSP 
can be obtained when non-universal third family sfermion masses are introduced. It 
should be noted, however, that none of these models combine a neutralino LSP with 
a long-lived stop as is done in the models studied here. 

There have been a number of searches [4] for i?-hadrons at colliders. The most 
recent searches were performed with data from LEP[17, 18] and the Tevatron[19, 20]. 
Guided by these results, we consider the possibility of detecting at the LHC stable 
stops with masses in excess of the excluded lower limit ~ 249 GeV/c^[21]. This work 
complements studies which have been made of the discovery potential of i?-hadrons 
at the LHC [22, 23, 24, 25] and at cosmic ray facilities such as the Pierre Auger [26] 
and IceCube[27] experiments. 



2. Long-lived stops in the MSSM-8 

A long-lived stop in models with a neutralino LSP can arise from a number of different 
parameter settings in the MSSM[28]. In such scenarios the stop is the NLSP and can 
be long-lived if it is very close in mass to the LSP. In that case only a very limited 
number of decay modes are open. In fact if the mass difference between the stop and 
the LSP is less than the IV-mass, only the suppressed modes: 

t^cx? (2.1) 
^ uxi (2.2) 
^ bffxl (2.3) 

are allowed. Figure 1 shows the Feynman graphs for the decay t — )■ ex?. 

In most phenomenological SUSY scenarios (such as MSSM-7[2] or mSUGRA/ 
CMSSM[30]), for simplicity one assumes a degeneracy in the sfermion sector (either 
at the Grand Unified Theory (GUT) scale, or at the electroweak scale). As we 
are here mostly interested in models with a light stop, but not necessarily with all 
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Figure 1: One loop stop decays [29]. 



sfermions being light, we shall break this degeneracy by introducing an additional 
mass parameter for the third generation of squarks. Our phenomenological MSSM-8 
model then has the following free parameters, 

• the top trilinear coupling At 

• the bottom trilinear coupling 

• the ratio of the vacuum expectation value of the two Higgs doublets tan P = ^ 

• the pseudoscalar Higgs mass 

• the fj, parameter 

• the common SU(2) gaugino mass M2 

• the first and second generation sfermion mass mo 

• the third generation sfermion mass mos 

This is the same model as the MSSM-7 model in Ref. [2], but with the extra param- 
eter which replaces rriQ in the (3,3)-component of the squark matriceses Mq, Mjj 
and Mr, (see Ref. [2] for more details). The non-universal sfermion mass scenario, 
introduced to separate the third generation of squarks from the first and the second, 
allows for the construction of models in which the stop and neutralino are close in 
mass. Note that a similar mass spectrum can be obtained in mSUGRA/CMSSM 
with a non-zero trilinear coupling at the GUT scale although here we prefer to work 
in the more general MSSM-8 model. Furthermore, we only need one gaugino mass 
parameter as we make the usual GUT assumptions for the gaugino mass parameters: 

Ml ~ O.5M2 (2.4) 

M2 ~ O.3M3. (2.5) 

The degeneracy between the stop and the neutralino masses influences the neu- 
tralino relic density via so-called coannihilations [31]. Hence, if we let the stop be 
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arbitrarily close in mass to the neutralino we would typically have too much coan- 
nihilation, leading to an effective annihilation cross section which is too large and a 
relic density which is too low. This condition is especially severe for lighter neutrali- 
nos, where the neutralino-neutralino annihilation cross section is fairly large already 
without coannihilations. Therefore, to respect the relic density constraint, we can- 
not allow mass differences between the stop and the neutralino which are too small 
at low masses, whereas at higher masses (TeV scale), we can let the stop be more 
degenerate with the neutralino. 

3. Model scan 

Models with a tiny mass difference between the stop and the LSP were generated 
using DarkSUSY[32]. The models were checked against current experimental con- 
straints, as implemented in DarkSUSY: 

• masses of superpartners [33] 

• masses of Higgs particles, checked with HiggsBounds 1.0.3 [34] 

• the invisible width of the Z boson [33] 

• the p-parameter [33] 

• the rare decay h — )■ S7, where we assume that the standard model and SUSY 
contributions gives a decay branching fraction in the range 2.71 x 10^"^ < 
Brih — )■ S7) < 4.39 x 10~^, in agreement with the experimental world average 
[35] and allowing for theoretical and experimental errors 

• the magnetic moment of the muon, where we assume that —2.7 x 10~^° < a^ = 
{g — 2)^/2 < 49 X 10^^*^, in agreement with the experimental measurement [36] 
and allowing for theoretical and experimental errors 

For more details on these limits, see Ref. [32]. 

In total, we have generated many thousands of different models with different 
masses and compositions. The focus has been on models with a small mass difference 
between the neutralino and the lightest stop. We have also chosen to apply the relic 
density constraint on Q^-h"^ with various degrees of hardness, 

1. a hard requirement 0.1099 ± 2 ■ 0.0062, corresponding to the five year WMAP 
result ±2(t[37]. 

2. a soft requirement 0.05 < l^x^^ < 0-2; corresponding to ^xh'^ = 0.1099l;]^oa- 

3. no requirement on the value of Q^h"^ ^■ 

^In practice all our models have a value of fl^h^ well below the five year WMAP result. 
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For practical reasons, we have selected a set of representative benchmark models 
with different neutralino masses by binning our models in 100 logarithmic intervals 
between 10 GeV/c^ and 10 TeV/c^ in neutralino mass. In each bin, we have selected 
the model passing all the constraints given above (for the different scenarios of the 
relic density) that has the smallest relative mass difference between the LSP and the 
lightest stop i.e. where (m^ — m^)/m^ is the lowest. This was done to ensure long 
stop lifetimes. 

The parameters and masses of the chosen models are written to SUSY Les 
Houches Accord 2 (SLHA)[38] files. These are used as inputs to SDECAY[39] 
which calculates the decay rate and width of the SUSY particles using loop correc- 
tions such as those shown in Figure 1. It is important for the decays in Eq. (2.1) 
and (2.2) that both the c and u masses are set to zero by default. The c-quark mass, 
rric = 1.21 GeV/c^, is therefore added to SDECAY by hand to avoid overestimating 
the stop decay rate contribution from Eq. (2.1). The w-quark mass is left as zero, 
which is a conservative choice; any variations of the mass would lead to a longer 
stop lifetime. SLHA output from SDECAY contains the decay branching ratios 
and widths for all sparticles. From the width, the stop lifetime can be obtained, 
and models that provide stops with a lifetime longer than 100 ns are then used as 
input to the event generator Pythia[3], version 6.4. The Pythia program generates 
events and calculates both the total MSSM cross section and the cross section for 
each sub-process. Kinematical event distributions such as momentum, angles and 
particle multiplicities are also provided. 

Many of our benchmark models depend sensitively on mass differences and hence 
on the input parameters. As the relic density depends exponentially on the mass 
difference between the lightest neutralino and the lightest stop in the coannihilation 
region, it is important to have sufficient numerical precision for this calculation. 
We have therefore generated all models and calculated accelerator constraints and 
the relic density in double precision. When we have calculated the decay rates 
with SDECAY, we have used SLHA2 to transfer our model parameters, which does 
not allow for double precision. Hence a small numerical error occurs. For a few 
models with extremely small mass differences, this numerical error was found to be 
significant, and those models were excluded from our set of models. For the remaining 
models, the numerical error is small enough not to be of any problem for the current 
study. 

Figure 2 (left) shows the lifetime of the stop versus the stop-neutralino mass 
difference and (right) the stop/neutralino mass difference versus the neutralino mass 
for the selected benchmark models. Figure 2(a) shows how the lifetime increases with 
decreasing stop-neutralino mass difference and also that for some models the mass 
difference is extremely small, at the order of 10~^ MeV/c^. The Figure also shows 
how the lifetimes range from around 10 s to 10^^^ s. We see that we can typically 
find models that pass all constraints for a variety of mass differences. Strong stop- 
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Figure 2: The left plot shows plots of the stop lifetime as a function of the mass 
difference between the stop and neutralino for all models. The right plot shows plots of 
the stop/neutralino mass difference versus the neutralino mass. 

neutralino mass degenracy on the other hand can only be found for neutralino masses 
above 200 GeV/c^. 

In more traditional SUSY models the stop is expected to be very short-lived 
with a lifetime around 10"^"^ s. From this perspective a majority of the models in 
Figure 2 give rise to long-lived stops. This study however focuses on stops which 
can both hadronise, traverse and leave a particle detector without decaying. The 
largest particle detector currently in use is ATLAS [23] which is 44 m long and has 
a diameter of 25 m. A particle moving with the speed of light will traverse the 22 
metres from the interaction point to the end of ATLAS in 73 ns. The stop masses 
in this study lies between 255 GeV/c^ to 2 TeV/c^ and can move with speeds less 
than the speed of light. A minimum stop lifetime is therefore chosen to be 100 ns, 
this lifetime allows particles with a speed v = 0.7c to traverse the maximum flight 
path in ATLAS. From a detector analysis perspective this is quite conservative since 
most particles will move in the plane transverse to the beam-line giving a flight path 
much shorter than 22 m. 

Figure 3 shows in the left column plots of the stop lifetime as a function of 
the stop-neutralino mass difference and in the right column plots of the neutralino 
mass as a function of the stop mass, after requiring the stop lifetime to be greater 
than 100 ns. The models are divided according to the hard (top), soft (middle) 
and no (bottom) ^^^^ cuts. The most striking tendency is how the number of 
models with a long-lived stop increases as the ^^h"^ cut is made looser. Also note 
from the plots 3(b) and 3(d) and 3(f) how loosening the cut on fi^^^ allows for 
lighter stop-masses. This is understood by recalling that the neutralino-neutralino 
annihilation cross section increases with decreasing neutralino mass and that the 
neutralino-stop coannihilation cross section increases when the stop and neutralino 
are close in mass. When respecting the relic density as measured by WMAP[37] the 
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Figure 3: The left column shows plots of the stop lifetime as a function of the mass 
difference between the stop and neutralino for the three different cuts on Q^h'^. The right 
column shows plots of the neutralino mass versus the stop mass. All plots show models 
where the stop lifetime is greater than 100 ns. 



contribution to the effective neutralino annihilation cross section from the neutralino- 
neutralino annihilation and the stop-neutralino coannihilation needs to be balanced. 
For this reason only TeV-scale neutralinos are allowed to be close to mass degenerate 
with the stops when ^^^"^ restricted, as a comparison between Figures 3(a) and 3(b) 
as well as 3(c) and 3(d) show. Note that this is also the case when we allow the relic 
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Figure 4: Inclusive stop pair-production cross section as a function of the stop mass, for 
the three Oj^/i^ requirements. 



density to be more than ±10a away from the central WMAP value, corresponding to 
the soft ^x^'^ requirement. When the relic density is unrestricted the contributions 
to the effective neutralino annihilation cross section also become unrestricted and 
low mass neutralinos can be close to degenerate with the stops. 

As seen from each of the plots 3(a), 3(c) and 3(e) the 100 ns cut leaves only 
models in which the neutralino and the stop are close to being mass-degenerate. The 
largest mass difference giving a long-lived stop is approximately 200 MeV/c^ and 
is found for the softer fi^^^ requirement in Figure 3(c). Some models with a mass 
difference at the order of MeV/c^ are also found in 3(a). However, most models 
with a long lived stop have a stop-neutralino mass difference below 0.1 MeV/c^. The 
extreme cases are found in Figure 3(e) where several models have a mass difference 
of the order 10~^MeV/c^. Note that no models where the decay t — ?■ ex? is open 
pass the 100 ns lifetime cut; all models passing this cut decay through i uxi- 

From this section it has become apparent that in order to get a long-lived stop 
in SUSY models with a neutralino LSP, the stop-neutralino mass difference has to 
be extremely small. Although it entails a certain degree of fine-tuning, approxi- 
mately mass-degenerate scenarios are chosen. The motivation for this is that all 
non-excluded regions in the parameter space need to be investigated and the chosen 
models would give rise to distinctive experimental signatures. 



4. Phenomenology of long-lived stop models at the LHC 

Figure 4 shows the inclusive stop pair-production cross section for the models with a 
stop lifetime greater than 100 ns. The calculation was made with Pythia. Relatively 
high cross sections of the order of 10 pb for pp-collisions at ^/s = 14 TeV can be 
achieved for low stop masses in the case where Q^h"^ is unconstrained, see Figure 
4(c). When Q^h'^ is restricted the maximum cross sections are 10~^ and 10~^ pb. 
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Figure 5: Higgs pair-production mechanisms. 



see Figure 4(b) and 4(a) respectively. In all three plots the cross section falls rapidly 
with mass and for masses above 1 TeV/c^ it falls below 1 fb. 

In models where the stop is the NLSP one would naively expect that pair- 
production of stops would dominate the Beyond the Standard Model (BSM) physics 
processes at a particle collider experiment. However, a closer look at the main 
production mechanisms at a 14 TeV proton proton collider shows that this is not 
necessarily the case. The SUSY models can be categorised in three cases depending 
on which is the dominant BSM process; 

1. Higgs production, mainly / + /'—)■ H^ + Hq and f + f^ H^ + Hq accompanied 
by the production of the other Higgs particles, see Figure 5. 

2. First and second generation squark production, mainly + qj qi,LR + Qj,LR, 
see Figure 6. 

3. Third generation squark anti-squark production, mainly g + g ^ ti + ti, see 
Figure 7 

Higgs-production: Figure 5 shows that the leading order Higgs pair production 
mechanisms happens via exchange oiWoiZ bosons in the s-channel. The lighter 
neutral Higgs mass m/jo is restricted from above and is expected be below 135 GeV/c^ 
in the MSSM [28]. The heavier neutral Higgs the pseudoscalar and the 
charged Higgs can be arbitrarily large. However, m^o and can be expressed 
as functions of m^o and will be strongly dependent on the value set for this parameter 
of the MSSM-8 model. All the heavier Higgs masses will in fact be very close to the 
value of m^o, usually not separated by more than a few GeV/c^. 
Production of squark-pairs: Figure 6 shows the squark pair production mech- 
anisms important in this work. Figure 6 describes the strong interaction processes 
and will in most cases dominate completely. Squarks could also be produced in elec- 
troweak processes through e.g. Higgs, gauge boson, neutralino or chargino exchange. 
However, the contribution from such processes is expected to be small compared to 
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Figure 6: Squark pair-production processes. 
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Figure 7: Quarlc-anti-quark anniliilation and gluon fusion diagrams producing a squark 
anti-squark pair. 



the strong production processes. We have verified that this is indeed the case for our 
models. 

Production of squark-antisquark pairs: Squark anti-squark pairs are produced 
via gluon-gluon fusion shown in Figures 7(b), 7(c), 7(e) and 7(f). Quark-anti-quark 
annihilation diagrams like the ones shown in Figures 7(a) and 7(d) also contribute. 
Even in this case, electroweak processes could contribute, but are expected to be 
subdominant compared to the strong processes. We have also here verified that this 
is the case for our models. 

In practice, we calculate the production cross sections with Pythia, including 
the strong processes given above [40]. 

Whether a scenario is dominant depends on the masses of the Higgs-particles, 
the squark sector and the gluino (recall that for all these models the LSP and the 
lightest stop is very close in mass). For models with relatively light squarks the strong 
processes in Figure 7 dominate. Stop-antistop pair production will dominate, since 
the first two generations squarks have considerably higher mass. For models where 
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Table 1: Masses of the stop, sup, pseudoscalar Higgs, and gluino as well as the main 
production mechanisms for models with ^^h"^ = 0.1099 ± 2 x 0.0062. 

the squark masses are large, the gluon-fusion process becomes suppressed by the low 
probability to obtain incoming gluons with a large fraction of the proton momentum. 
Instead, for intermediate values of the squark masses, the strongly coupled process 
of Figure 6 becomes leading. In this case the flavour of the proton valence quarks 
translates into production of squark-pairs from the first two generations. For models 
where either the squark-masses or the mass of the gluino-propagator gets very large, 
weak production of Higgs-bosons, as in Figure 5, will dominate. 

Table 1 lists the seven models with a long-lived stop and Q^^.h'^ = 0.1099 ± 2 ■ 
0.0062. All models are given a bookkeeping number and the values of the stop'^, 
sup(Mi:,), pseudoscalar Higgs and gluino mass are listed. The models are divided 
into three cases according to which BSM production mechanism dominates. Case 1 
represents Higgs production, case 2 first and second generation squark production 
and case 3 third generation squark production. For each main group of production 
mechanisms the four largest contributing processes are listed in the last column of 
the table. 

As seen from Table 1 the BSM contribution for five of the models are dominated 
by Higgs production, one produces mainly squarks and one stop-anti-stop pairs. The 
five models which are dominated by Higgs production also have a large mass dif- 

^When referring to the stop, we here mean the lightest stop, which we denote ti. 
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Table 2: Masses of the stop, sup, pseudoscalar Higgs, and gluino as well as the main 
production mechanisms for models with 0.05 < Q^h^ < 0.2. 

ference between the squark and the Higgs sector. Also note that the gluino mass 
for all five models is above 10 TeV/c^, which strongly suppresses the strong produc- 
tion of squarks that otherwise would dominate. This explains why the electroweak 
production of Higgs particles dominate the overall BSM particle production. 

Table 2 shows the main production mechanisms for long-lived stop models with 
0.05 < f^x/i^ < 0.2. Higgs production dominates six models, four are dominated by 
first and second generation squark production and three models are dominated by 
third generation squark production. The underlying mechanisms determining which 
types of process dominate are the same as before. 

When fij^/i^ is unrestricted production of stop and sbottom dominates, closely 
followed by production of first and second generation squarks and a small fraction 
of models dominated by Higgs production, see Table 3. The tendencies discussed 
before become even clearer for this choice of VL^h"^. Higgs and squark production 
dominate when the stop mass is quite high. Higgs production is prominent when 
Wjao is relatively small and squark production is suppressed due to heavy masses. 
When the Higgs masses increase and the squark masses decrease and the first and 
second generation squarks are comparable in mass to the third generation, the model 
is dominated by squark production. For models with smaller squark masses there 
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10 


782 


883 


983 


5409 


1i ~ 




' qi,L 1 






11 


703 


767 


1030 


4865 


y 






— n ■ T 




12 


890 


1020 


889 


6160 












13 


430 


464 


1183 


2976 












14 


346 


4754 


603 


2392 












15 


542 


2523 


352 


3749 












16 


498 


2728 


1321 


3447 












17 


463 


907 


1162 


3203 












18 


1073 


4483 


1181 


7424 


g- 


^g 


^ ti 4 


-h 


3 


19 


604 


974 


885 


4176 


f- 




^ ii 4 


-h 


20 


756 


2885 


1026 


5231 


g- 


^g 




-k 




21 


1110 


1528 


1107 


6988 


g- 


Kg 


^2 ^ 






22 


653 


1242 


363 


4516 












23 


255 


4655 


1343 


1766 












24 


288 


3030 


1268 


1993 












25 


402 


2910 


1103 


2782 











Table 3: Masses of the stop, sup, pseudoscalar Higgs, and gluino as well as the main 
production mechanisms for models with no restrictions on O^/i^. 

can also be stop-anti-stop production through gluon-gluon fusion. Third generation 
squark-anti-squark production dominates the BSM production when the stop is much 
hghter than the other squarks and the heavy Higgs particles. However, as long as 
it is light enough for the gluon-gluon fusion processes to be the most important 
production mechanism it will still dominate the overall BSM processes. 

From Tables 1-3 we can deduce that models containing a long-lived stop with 
mass above 1500 GeV/c^ could be very hard to detect at the LHC Not only because 
the cross section falls rapidly with mass, but also because these models are dominated 
by Higgs production. In fact for some of these models it will be difficult to prove 
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that supersymmetric processes are at work since none of the standard SUSY candles 
hke high leptons or missing transverse energy will be present. The Higgs particles 
produced will decay to bb or ti and if the first decay dominates it will be hard to 
distinguish from background. For some of these models it could also happen that 
a long-lived stop signature will be the first indication of the presence of beyond the 
Standard Model physics. 

Models dominated by first and second generation squark production can also be 
quite challenging since the squarks will rapidly decay through the mode qi Xi + Qi ■ 
These decays will give rise to hard jets and missing transverse energy. This signature 
will surely deviate from standard model processes, but due to the high squark masses 
in these models discovery will require 14 TeV collisions and high statistics. In most of 
these models approximately 6% of the total BSM processes produced in pp collisions 
are stop-pair-production. For this reason and of course depending on the mass and 
the cross section there should be a possibility of discovering such states. 

Long-lived stops are of course easiest to discover in models where stop production 
dominates the overall SUSY processes. The limitations when it comes to discovering 
such states arise from the cross section and how well the signal signature can be 
separated from the Standard Model background. This will be studied in more detail 
in the following sections. 

5. Discovering long-lived stops at the LHC 

If long-lived stops exist it will be of great importance to discover them at the LHC. In 
the following sections we will use the performance of the ATLAS [23] detector as an 
example of a general purpose detector and investigate the possibility of discovering 
long-lived stops using generic and global cut variables. All events are generated using 
Pythia 6.4. 

5.1 Interactions of long-lived stops 

Long-lived stops produced at a hadron collider will immediately form i?-hadrons. 
Both i?-mesons tiq and R-baryons tiqq can be formed. We adapt the assumption 
that 90% of the stops will initially form i?-mesons and 10% form i?-baryons from 
Ref. [41]. Through interactions with the detector material most mesons convert to 
baryons. In a conservative assumption any of the baryon states tiuu (charge= +2e), 
iiud {+e) or iidd (Oe) are equally likely for a stop exiting the detector. The situation 
for the anti-stop is somewhat different; it will also form an i?- meson when produced, 
but a transition to a anti-baryon state like tiuu is unlikely and if it occurs it will lead 
to an immediate annihilation back to a meson. An anti-stop exiting the detector 
will therefore be in a mesonic tiu or iid state, with charges and — e, respectively. 
Given these assumptions, roughly 80% of stop events are expected to have at least 
one charged i?-hadron leaving the detector. 
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Model 


GeV/c^ 


(j~ = 
(fb) 


N(100 fb-i) 




1 


255 


15000 


1.5-106 


0.0019 


2 


542 


306.6 


30660 


0.0047 


3 


746 


39 


3900 


0.0128 


4 


921 


9.7 


973 


0.0614 


5 


1073 


3.9 


390 


0.0442 


6 


1588 


0.22 


22 


0.1134 



Table 4: Inclusive stop pair cross section, expected number of events at 100 fb and 
O^^./i^ for the six selected SUSY models. 

Due to ionisation energy losses, the charged i?-hadrons will leave tracks in the 
detector's tracking devices. Very little energy will be deposited in the calorimeters 
because of an expected small interaction cross section[41, 42, 43, 44]. The interactions 
that do take place here are responsible for the meson to baryon conversion which can 
also lead to a change of electrically neutral states to charged states and vice versa. 
Finally, charged i?-hadron states also leave tracks in the outermost parts of a detector 
which are designed to detect muons. 

An i?-hadron remaining charged when traversing a detector will have a similar 
signature as a high transverse momentum (pr), possibly slow moving, muon. Jet 
activity in the event is expected to originate from initial/final state radiation or from 
stop production through decay of heavier sparticles. Since there is a probability for 
the i?-hadron to acquire or lose its charge, there can also be events where a high 
inner detector track disappears before reaching the muon detectors, or that a high pt 
muon track has no inner detector track matching it. These effects are best studied 
with more detailed detector simulations as in Ref . [23] , and are outside the scope of 
this paper. 

5.2 Selection of benchmark models 

A set of six models with stop masses ranging from 255 to 1588 GeV/c^ was selected. 
The mass spectra and some details of these models are listed in Tables 4 and 5 
respectively. 

The models were selected to provide a continuous and wide mass range. All 
samples are scaled to 100 fb~^ corresponding to one year of running 14 TeV collisions 
at a luminosity of 10~'^'^cm~^s~^. Some models, for example those containing the 255 
and 542 GeV/c^ stops, have, in principle, large enough cross sections to make earlier 
discoveries possible, although this has not been investigated in this work. 

All selected models have third generation squark production as the main BSM 
production mechanism. The total fraction of stops produced in a given model de- 
pends not only on the direct production of stop pairs but also on the indirect pro- 
duction via decay of heavier sparticles. The largest contributors to the latter are the 
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Model 


tl 

GeV/c2 


Percentage of total / Fraction decaying to titi 




hh 


hh 


^262 


1 


255 


92% 


1 


0.3% 


0.36 


6 % 


0.48 


1 % 


0.77 


2 


542 


72% 


1 


2% 


0.36 


13% 


0.96 


10% 


0.98 


3 


746 


39% 


1 


11% 


0.16 


35% 


0.00 


14% 


0.74 


4 


921 


23% 


1 


15% 


0.00 


23% 


0.32 


20% 


0.00 


5 


1073 


38% 


1 


12% 


0.19 


34% 


0.00 


15% 


0.81 


6 


1588 


30% 


1 


11% 


0.62 


20% 


0.00 


18% 


0.92 



Table 5: The percentage of third generation squark production contributing to the total 
BSM cross section combined with the fraction of stop and sbottom events decaying to titi. 



Process 


a 
(pb) 


N(100 fb-i) 
10^ 


Z+ jet 


138.4 


14 


t^+jet 


1192 


119 


WW, ZW, ZZ 


231 


23 


tt 


490 


49 


bb 


5063 


506 


QCD J5 


12490 


1249 


QCD J6 


359 


36 


QCD J7 


16 


2 


QCD J8 


0.024 


0.0024 



Label 


Pt interval 

(GeV/c) 


QCD J4 


140-180 


QCD J5 


280-560 


QCD J6 


560-1120 


QCD J7 


1120-2240 


QCD J8 


> 2240 



Table 6: Background cross sections and 
expected number of events at 100 fb~^. 



Table 7: pf-range of the hard scatter- 
ing process defining the QCD kinematic 
slices. 



production and decays of the other third generation squarks. These are produced 
via the processes g + g ^ Qi + Qi and f + f — )■ Qi + Qi, where Qi denotes ii, ^2, 
bi or 62- Table 5 shows the contribution from each of the third generation squark 
production processes to the total BSM cross section. The fractional numbers repre- 
sent, in the cases of t2 and 61^2, the decays to a tiii pair. For example, in model 2 
bibi production will constitute 13% of the total BSM cross section and 96% of all 
bibi pairs will decay to tiii. Based on these numbers one can calculate the inclusive 
iiti cross section (c^jj and then get the expected number of iiti pairs at 100 fb~^ 
(N(100 fb~^)) as shown in Table 4. The last column of Table 4 list the value of ^^^"^ 
for each model, which in all cases except for model 6 is well below the limit computed 
from the WMAP- measurement. 

5.3 Standard Model background processes 

All Standard Model processes giving rise to high-p^^ muons could mimic the signal 
of a long-lived, hadronising stop. The most important processes are: 

a) Z + jet 

b) W + jet 
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c) WW, ZW, ZZ (diboson production) 

d) ti 

e) bb 

f) high pt QCD-processes not covered in the above hst. 

Muons are produced in the processes Z — )• fi^fi~ and W"^ — )■ jju^u^ip^). In 
the 2—7-1 processes producing single vector bosons, the muons typically carry a 
momentum which is half of the boson mass. These muons do not have high enough 
Pt to constitute any real background for the long-lived stops. In the 2 — )■ 2 processes 
where the boson is recoiling against a jet, the boson itself carry momentum which 
is transferred to its decay products. This allows for high-p^^ muons coming from Z 
and W decay and is the reason for the choice of a) and b). 

The top quark decays to a.W and 6-quark and gives rise to muons either through 
the decay of the W or the decay of a S-meson, or both. In bottom pair-production 
all muons arise from the decay of 5-mesons. 

Background f) is referred to as high-p^ QCD-processes. These processes produce 
the light quark flavours and/or gluons giving rise to multiple jets. Muons arise 
primarily from decays of charged pions and kaons. These processes will be dominating 
at the LHC and in order to generate the number of events needed to ensure reasonably 
small statistical uncertainties it is common to divide events into pt intervals according 
to the hard scattering. Here we use the same division as in Ref. [23] and the intervals 
are listed in Table 7. Only the samples J5-J8 will be used in the following analysis 
as it turns out that the J4-background is negligible. 

5.4 Data selection and detector acceptance cuts 

A candidate event sample was created by requiring at least one i?-hadron / muon with 
transverse momentum larger than 100 GeV/c. The pr spectrum of the i?-hadrons 
and muons after this preselection is shown in Figures 8(a) and 8(b) respectively. 
The preselection removes a substantial fraction of the background, while most of the 
signal is kept, which is also seen in Table 8. As expected the signal retention is larger 
for high masses. The exception for the 921 GeV/c^ stop mass case is explained by 
Table 5 where we see that 35% of the total SUSY production give rise to events 
where no stop pairs are produced. Out of the events in this sample that contain 
stops 98% pass the pre-selection. 

i?-hadrons/muons are required to be within the pseudorapidity range |?7| < 2.4. 
This requirement is motivated by the geometrical acceptance of the ATLAS first 
level muon trigger system, which extends from 77 = to 77 = ±2.4. The signal and 
background distributions of t] is shown in Figures 8(c) and 8(d). 

Since i?-hadrons are quite heavy they can move with a speed /3 = ^ < 1. As the 
particles need to reach the trigger stations within 25 ns to be associated with the 
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Figure 8: The transverse momentum, r] and f3 distributions for signal i2-hadrons 8(a), 
8(c), 8(e) and the various muon backgrounds 8(b), 8(d), 8(f). All plots show the spectrum 
after a preselection including only events with at least one i?-hadron/muon that has a 
transverse momentum larger than 100 GeV/c. 



correct bunch crossing, the requirement of /3 > 0.7 is applied. The /^-distributions for 
signal and background are shown in Figure 8(c) and 8(f) respectively. Comparison 
of the two figures shows that the /3-spectrum for the signal tends to peak at lower 
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Figure 9: Summed transverse momentum of tracks in a cone AR < 0.2 and the distance 
AR between i?-hadron/muon track and the closest jet for signal -R-hadrons 9(a), 9(c) and 
the various muon backgrounds 9(b), 9(d). All plots show the spectrum after a preselection 
including only events with at least one i?-hadron/muon that has a transverse momentum 
larger than 100 GeV/c. 



values for higher masses, while the background sample has nearly all events at /3 = 1. 
This is confirmed by Table 8 showing that 100% of events in each background sample 
pass this cut and that the fraction of signal events passing the cut decreases rapidly 
with mass. 

In addition to the pre-selection and detector acceptance cuts, particle isolation 
cuts are imposed. Firstly we require that the sum of the transverse momenta of 
tracks within a distance AR = (Ary)^ + (A0)2 < 0.2 of the R-hadron/muon track 
is less than 10 GeV/c. Figures 9(a) and 9(b) show the summed of tracks within 
a cone of size 0.2 around the i?-hadron/muon track. The background distribution 
(Figure 9(b)) shows that muons are very often produced close to other high pj- tracks, 
while this is to a much lesser extent the case for the signal, ~90% or more of the 
QCD background is removed by this cut, cf. Table 8. Most of the muons remaining 



-19- 



Sample 


Pre- 
selection 


\r]\ < 2.4 


/? > 0.7 


/ J A R^n 9 Pt 

< lOGeV/c 


AR > 0.4 


Combined 


255 


77 


97 


85 


100 


96 


73 


542 


93 


98 


73 


100 


96 


61 


756 


97 


99 


65 


100 


98 


54 


921 


64 


99 


61 


100 


99 


52 


1073 


98 


99 


56 


100 


98 


47 


1588 


97 


99 


43 


100 


98 


37 


QCD J5 


0.17 


99 


100 


12 


97 


3 


QCD J6 


0.98 


100 


100 


9 


95 


2 


QCD J7 


0.97 


100 


100 


7 


93 


2 


QCD J8 


1 


100 


100 


5 


90 


2 


tt 


2 


98 


100 


93 


67 


58 


bb 


0.08 


97 


100 


5 


99 


1 


Z + jet 


15 


89 


100 


100 


95 


82 


W + jet 


4 


83 


100 


98 


90 


72 


WW,WZ,ZZ 


0.14 


85 


100 


99 


96 


79 



Table 8: Percentage of remaining events after a given cut. The pre-selection column is the 
fraction of events remaining after requiring at least one muon with transverse momentum 
larger than 100 GeV/c. The remaining five columns show the efficiency (in percentage 
points) of each cut separately, relative to the preselected sample. 

after these cuts originate from the decay of W ot Z bosons. 

The other isolation criteria requires that the distance AR between the muon and 
the closest jet is larger than 0.4. This result is highly dependent on the jet-algorithm 
used; in this case jets are defined using a cone algorithm with cone size = 0.7 and a 50 
GeV energy threshold. The distance between a i?-hadron/muon and the closest jet 
is plotted in Figures 9(c) and 9(d) respectively. Both distributions show a relatively 
good separation between i?-hadrons/muons and jets, which is also reflected in the 
second last column of Table 8 which shows that more than 90% of the events are 
kept for both signal and background. 

The last column of Table 8 shows the total remaining percentage of the signal 
and background samples after applying all cuts. The cuts remove between 30% and 
60% of the signal. The requirement of /3 > 0.7 is the one which removes, by far, the 
most signal. Out of the background events up to 99% of the QCD and bb samples are 
rejected by the combined cuts and it is the muon isolation criterium that removes 
most events. In the case of electroweak processes ( ^+jet, H^+jet and WW,ZW,ZZ 
) the samples are predominantly reduced by the cut on |?7| and the distance between 
the muon and the closest jet. The latter requirement also removes the bulk of the 
rejected tt events. 

5.5 Enhancing the stop i?-hadron signal 

After the acceptance cuts and the isolation requirements on the muon candidates we 
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Figure 10: Leading jet transverse momentum distribution for signal 10(a) and background 
10(b). 



Sample 


Number of Events 


Pre- 


Isolation and 


Jet pt 


Jet and /i 




selection 


Acceptance 


< 100 GeV/c 


Pt -cut 


255 


1.02885e+06 


753560 


485212 


18234 


542 


24019 


14578 


7702 


2582 


756 


5946 


3223 


1860 


1083 


921 


1890 


978 


643 


438 


1073 


582 


272 


156 


127 


1588 


33 


12 


6 


5 


j5 


2.17032e+06 


62999 


4734 





j6 


353081 


8100 


127 





j7 


5556 


97 


1 





j8 


3 


0.04 








tt 


1.07456e+06 


626064 


105812 


24 


bb 


450008 


5971 


910 





Z + jet 


2.03223e+06 


1.66846e+06 


518456 


99 


W + jet 


4.68819e+06 


3.35225e+06 


992362 


34 


WW,WZ,ZZ 


32469 


25701 


12586 


18 



Table 9: Number of events remaining after preselection cut, isolation and acceptance 
cuts, leading jet cut and leading jet + leading muon cut. 

are left wfth from 12 to order of 750 000 events for the signal samples, while the total 
number of background events is over 5 million events. The W/Z + jet events are the 
dominant background followed by QCD and tt events, as shown in Table 9. 

After the preselection requirement of a transverse momentum on the muon can- 
didate [prifJ') > lOOGeV/c) the main contribution from the electroweak processes 
comes from bosons produced in 2 — t- 2 processes in which a jet recoils against the 
boson. This gives a boost to the boson which leads to high px of the muons produced 
when it decays. 
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Figure 11: Transverse momentum distribution for all i?-hadron signal samples and the 
summed muon background after applying preselection, detector acceptance and isolation 
cuts and requiring that pT < lOOGeV/c for the leading jet. 



In the signal samples jets can be produced in various ways, in the case of direct 
stop pair-production jets originate from initial and final state radiation as well as 
the underlying events. When stop pairs are indirectly produced through the decay 
of heavier squarks, jets originate from hadronic decays of W or Z bosons produced 
in, for instance, sbottom decay (6j — )■ ti + W) or stop decay (^2 —> ii + Z). 

In the models considered here however, the mass-splitting between the third 
generation squarks is relatively small, leading to rather low momentum for the ac- 
companying jets in these decay chains. The distribution of the momentum of the 
leading jet in the signal and background samples, shown in Figures 10(a) and 10(b) 
demonstrates the expected difference between signal and background. 

To exploit this difference to further enhance the signal we require the leading jet 
to have transverse momentum below 100 GeV/c which leaves the number of events 
shown in column 4 of Table 9. The W and Z samples are reduced by 70%, while the 
corresponding numbers for the signal vary between 40% and 50 %. Also note that 
both the tt and bb samples as well as the QCD samples are significantly reduced. 

Figure 11 shows the transverse momentum distribution for all signal samples 
and the summed muon background after applying the cut on the px of the leading 
jet. The remaining background decreases rapidly from approximately 100 GeV/c 
and motivates an additional hard px-cut. 

To further reduce the remaining background a cut on the leading i?-hadron/muon 
is imposed. We plot S/ \/B + S as a function of the cut on muon-p^ to find the 
cut value that optimises the signal to background ratio, see Figure 12. Based on 
these plots a requirement of pt > 500 GeV for the leading muon is imposed. The 
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Figure 12: S/^JB + S for all signal samples as a function of the cut on the leading muon 
transverse momentum. 



remaining events are listed in the last column of Table 9 and from these numbers we 
can easily calculate the signal significance S/ \/ S + B which is larger than 5 for all 
signal samples except the heaviest mass point at 1588 GeV/c^. 
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6. Conclusion 



In this work we have studied the possibihty for a long-hved stop to arise in super- 
symmetric models with a neutrahno LSP. This is achieved when the mass difference 
between the stop and neutrahno is extremely small, less than the c-quark mass. The 
small mass difference increases the coannihilation cross section leaving the neutrahno 
relic density smaller than the WMAP-result. This implies either that the neutrali- 
nos coming from these models can only constitute a part of the total dark matter 
and that there in addition must be other particles/mechanisms making up rest, or 
that the lower bound on the relic density, as derived from the WMAP data, is in- 
valid due to modifications to the expansion rate of the universe prior to Big Bang 
nucleosynthesis as proposed e.g. in Ref. [45]. 

We have also seen that SUSY models with a long-lived stop give rise to non- 
traditional SUSY signatures. In many of the studied models the main processes will 
be various kinds of Higgs production or stop-pair-production, none of which give rise 
to the standard leptons plus jets plus missing transverse energy signatures. If such 
models are realised in nature the discovery of a long-lived stop could be the first sign 
of physics beyond the standard model. 

Finally we have studied six SUSY models where stop pair-production dominates. 
By using a range of final state observables, it should be possible to discover stable 
stops with masses up to around 1 TeV/c^ at the LHC. Should timing-specific observ- 
ables, such as the reconstructed speed of a i?-hadron candidate, be used in a search, 
the mass reach could be higher still. 
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